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Abstract 



D \ We analyze in detail the renormalization group flows which follow from the 

recently proposed all orders (3 functions for the Chalker-Coddington network 
model. The flows in the physical regime reach a true singularity after a finite 
scale transformation. Other flows are regular and we identify the asymptotic 
directions. One direction is in the same universality class as the disordered 
^ \ XY model. The all orders P function is computed for the network model of the 

O \ spin Quantum Hall transition and the flows are shown to have similar prop- 

erties. It is argued that fixed points of general current-current interactions 
.^ ■ in 2d should correspond to solutions of the Virasoro master equation. Based 

^ , on this we identify two coset conformal field theories osp{2N\2N)i/u{1)q and 

0^ I osp{4N\4N)i/ su{2)o as possible fixed points and study the resulting mul- 

^^ I tifractal properties. We also obtain a scaling relation between the typical 

t^^ ■ amplitude exponent qq and the typical point contact conductance exponent 

Xf which is expected to hold when the density of states is constant. 
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I. INTRODUCTION 

The Chalker-Coddington network model is perhaps the simplest model believed to be in 
the universality class of the Quantum Hall transition |[l[] . The model can be mapped to a free 
Dirac fermion theory with three independent random potentials [0 . The renormalization 
group (RG) fixed point of the disorder averaged effective action is expected to capture the 
critical properties of the transition. 

Recently, an all orders (3 function was proposed for the network model |Q. The main 
difficulty encountered in interpreting the resulting RG flows is the occurance of singularities 
in the (3 function. For the case of anisotropic sm(2), by using duality and certain topological 
identifications it was shown in ||^ that all RG trajectories were regular in the sense that all 
flows could be extended to arbitrarily large and small length scales without encountering any 
true singularities. Further checks of the (3 function in the case of anisotropic su{2) described 
in [Q provided additional support that the (3 function could be analytically continued beyond 
the domain of convergence of the perturbation series. 

The main purpose of this work is to analyze in detail the RG flows of the Chalker- 
Coddington network model and also the network model for the spin Quantum Hall transition 
P-H] and to attempt to resolve the possible singular flows using the ideas in |^. There are 
large classes of flows that are regular, and for these we identify two universal attractive 
directions at large distances, referred to as dXY and gX below. (Some of the phases listed 
as distinct in ^ become identifled since we now understand how to cross the singularities.) 
The dXY phase is in the same universality class as a disordered XY model or Gade-Wegner 
class 1^ and is equivalent to the model studied in ^^. Unfortunately the flows in the 
physical domain of the network model with all couplings (which are variances of disordered 
potentials) positive are pathological in the sense that they encounter a true singularity after 
a flnite scale transformation and the flow cannot be continued to larger length scales. These 
singular flows seem to be related to the fact that there are level current algebras in these 
theories which arise naturally in the supersymmetric method for disorder averaging. 

Putting aside the /? function, we argue that the possible flxed points of general anisotropic 
left/right current- current interactions should correspond to the solutions of the Virasoro 



master equation |TT| which are built upon the same tensors that deflne the current in- 
teractions and thus preserve the global symmetries. We classify these solutions for both 
network models and show that they correspond to current- algebra cosets. Since this argu- 
ment is independent of the (3 function, it remains a possibility that a proper resolution of 
the singular flows could lead to these flxed points. It was argued in [|] that such cosets can 
arise as flxed points of the RG when the coupling to a current subalgebra flows to infinity 
and the current subalgebra is gapped out. Some of the regular flows actually realize these 
flxed points, in particular phase C in the unphysical regime of spin network model flows to 
osp(4|4)i/osp(2|2)_2. By doubling the degrees of freedom as in |]TD|], this may correspond to 



a physical model of interest. For the spin network model one solution corresponds to the 
coset osp(4A^|4A^)i/s-u(2)o, where N is the number of copies, and this flxed point gives some 
of the same exponents as percolation [Q though the equivalence of the two theories is still a 
matter of debate. 

For the m(1) network model the only zeros of the (3 function correspond to a purely 
random u{l) gauge fleld. We point out that the stress tensor of this line of flxed points 



corresponding to the strength of the disorder does not satisfy the Virasoro algebra, which 
apparently has not been noticed before. There is a special point on this line that leads to a 
constant density of states, and can formally be viewed as "osp(2A^|2iV)i/M(l)o" as far as the 
computation of critical exponents is concerned. (This is not a true coset since the u{l)o does 
not commute with the stress tensor, however we find the coset notation convenient and will 
use it in the sequel.) This c = conformal field theory is different from the c = —2 theory 
osp{2N\2N)i/u{l) ^u{l) considered in |T^, nor is it the same as the theories considered in 



ME 

We wish to emphasize that due to the singularities of the RG flow in the physical regime 
of the network models we cannot show that the above cosets genuinely arise as fixed points. 
However, based on our analysis of the Virasoro master equation, they remain as interesting 
candidates for the correct fixed point and are worth further investigation, even with the 
aim of ruling them out. In this paper we thus also explore the multifractal properties of the 
above coset theories, and show that these cosets lead to multifractality in a very natural way. 
This multifractality is an important constraint on the critical theory and since it is not a 
generic feature of the kinds of fixed points that can arise, it can serve as a criterion for ruling 
out some theories. In particular we show that the coset osp{2N\2N)i/u{l) ® m(1) does not 
exhibit multifractality, whereas osp{2N\2N)i/u{l)o does. For the coset osp{4:N\4:N)i/ su{2)o 
we show that the typical amplitude exponent a^ = 9/4, and the typical point contact 
conductance exponent is X^ = 1. Since the map to percolation [§ does not apply to 
averages of higher moments of correlation functions, this result can serve as a useful test of 
the validity of the fixed point osp{4:N\4N)i/ su{2)o. For the u{l)o coset, it turns out ao = 3. 
However if one associates ao with the twist fields which modify the boundary conditions 



then oq is again 9/4 and this is close to numerical estimates of 2.26 |T^,[T6|. The small errors 



in the most recent simulations however seem to rule out 9/4 |T^. We also obtain a scaling 
relation between a^ and the typical point-contact conductance exponent Xt which we expect 
to hold in a theory with a constant density of states: Xt = 2{ao — 2) |T^. 

Our results are presented as follows. In section II we extend the computation in 0] to 
A^ copies, and show, as expected, that the P functions are independent of A^. We present 
a remarkable strong-weak coupling duality of the (3 function. We then describe all possible 
flows, and describe the nature of the singular flows. In section HI we extend this analysis 
to the network model for the spin Quantum Hall transition and find similar properties. In 
section IV we present our results on the multifractality. The subsections on the Virasoro 
master equation and multifractality are logically independent of the sections describing the 
RG fiows and can be read separately. 

II. THE CHALKER-CODDINGTON NETWORK MODEL 

The Chalker-Coddington network model [||] is expected to be in the same universality 
class as the Dirac Hamiltonian with random mass, potential and abelian gauge fields [Q : 

H- f ^ + '^ -idj + A-\ 

\-td, + A, V-m J ^^-^^ 

where A, m, V are centered gaussian random potentials with variance ga, Qm and g^ respec- 
tively. The physical regime of the network model thus corresponds to all g positive. The 



couplings ga,gv,gm correspond to randomness in the link phases, the flux per plaquet and 
the tunneling at the nodes respectively. This model also arose in the work ||^. The case 



dm = dv = was studied in 0,0,|I^, and the case Qm + Qv = in |]T0[. The one- loop 
j3 function for all y's non-zero was computed in [|l9l. The all-orders (3 function was proposed 
in§. 

A. N-Copies effective action. 

Following the conventions in |l^3 we study the action 



27r 



ij_{d,-iA,)ij^ + ij4ch-iA^)ij+ (2.2) 



where z = x + iy,z = x — iy. The above action corresponds to S* = i Jip^Hip where H is 
the hamiltonian of the network model. To study the Green functions at energy S, one lets 
H —>■ H — S which leads to a term in the action 

Se = if-^£^£, $^ = ?_^+ + 7^-?+ (2.3) 

J Zn 

The retarded and advanced Green functions then correspond to S = E±i£ with e small and 
positive. 

To study the delocalization transition one is interested in the critical points of the disorder 
averaged theory at £^ = 0. Since $£■ corresponds to various mass terms, the theory can not 
be critical when S ^ 0. The (3 function we determine below is for S = 0, and for this 
reason it is independent of whether the various copies are retarded or advanced. One can 
in fact convert advanced copies into retarded ones by flipping the sign of the right-moving 
fields ipj. -^ — ^±- Since this does not change the ip.^ operator products, the /3 function 
is unchanged. Once one flows to a critical point, one needs to reintroduce S in order for 
instance to compute conductivities. This can be compared with the sigma model approach. 
The latter approach was applied to the coupling sub-manifold g^ + gv = ^ oi the network 
model in [|lOl. An important point is that the saddle point approximation that leads to a 
sigma model occurs in the S = theory. Thus one expects that our RG flows should in 
principle capture these saddle points. 

In order to study disorder averages of mixed products of retarded and advanced Green 
functions we introduce N copies of the action ( |2.2|) , with fields ip^, ^± where a = 1, ..,N 
is a copy, or flavor, index. The disorder averaging is then performed using the super- 
symmetric method by introducing ghosts (3'^, (3^. To present the effective action in a 
compact form, we introduce the 2N component vectors \l/^, a = 1,..,2N, and we let 
\&± = ('j/'^, ...,'?/'^,/5^, ...,/3^) and similarly for \1/-|-. We define a grade [a] =0,1 where 
[a] = 1 for the fermionic components ip and [a] = for the bosonic components (3. The 
effective (disorder averaged) action is then 

S = Si,ee + j ^ {9+0+ + g.O- + gaO'') (2.4) 



where g± = g^ ± Qm-, and S'free the free c = conformal field theory of the ip^ j3 fields, 

5free = / ^ (^-5^^+ + ^-5.^+) (2-5) 

where \1/_\1/+ = X]a^-^+5 ^tc. The operator product expansion (OPE) imphed by the free 
action ( p.5| ) is 

^+(;2)^i(0) ~ -<5"^ *^(^)^+(0) ~ -tll^^ft^ (2.6) 

and similarly for ^ . In order to compute the (3 function we express the perturbing operators 



in eq. (2^) as left-right current- current interactions. Define the left-moving currents 

and similarly for the right movers J ^ and H . The perturbing operators can then be 
written as 

1 1 2iV 

0+ = '- {(^.^^f + (vi/_^+)2) = i 5: (-)W (jf j^_" + rj>i':) 

2N 

Q- = (^_^+)(^_^+) = Y. (-)["! //"^^ff^'* (2.8) 

a, 6=1 

2iV 
Qa ^ (^_^^)(^_^^) = ^ ^_)[a] + [b]H''aJjbb 

a,b=l 

The currents ( |2.7| ) satisfy the osp{2N\2N)k current algebra with level k = 1. 

B. /3eta functions and duality. 

We now extend the computation in [Q to A^-copies by using the previous compact no- 
tation. The resulting (3 function is independent of A^. We then describe the strong-weak 
coupling duality of the (3 function that is important for extending the flows to all scales. 

The effective action can be expressed as 

5 = 5^L + / ^ E 9aO^ with O^ = d^rl' (2.9) 

where S^'^^ is the action for the conformal WZW model with current algebra symmetry Qk, 
where k is the level, and J" are the Qk currents. For the network model Qk = osp{2N\2N)i, 
and the tensors d^^ are implicitly defined by eqs.(|2.8|). The /5 function proposed in [^ is 



expressed in terms of some OPE coefficients C, D, C. Let T"^ be the left-moving operator, 

T\z) = di,r{z)j\z) (2.10) 

Then the RG data can be computed from the OPE's 



C^(z,z)O^(0) ~ 4Cc'' C^(0) (2.11) 



T^(z)O^(0) ~ \ {2kD^'' + C^^) 0^(0) 



Specializing to the A^-copy network model one finds that C,D,C are independent of 
N due to the fact that Y.a=ii~y"'^ = O5 which is equivalent to the statement that the 
superdimension of osp(2A'^|2A^) is zero. As expected, the (3 function is then identical to the 
A^ = 1 result in [il: 



A 



8g+{2g,{gl + A) + {2-g_){gl + 2g_ 



g+ 



{4-gl){2-g,y 
8 gl{2 + g ^2 

(4-^ 
4 ({gl - g^_)il6 - glgl) + 4^,(?^(2 + g,){2 - g^f 



/3.- = T; "V (2-12) 



^9. 



{A-glf{2-g_y 



Here (3 = dg/dr where r = logr is the RG 'time' and r is a length scale. The only zero of 
the (3 function is at 5^^ = with ga arbitrary. We remark that (3g_ is always positive for real 
couplings so that g^ is always increased by RG transformations. 

The strong-weak coupling duality of the (3 function is of the following kind. Let g* denote 
dual couplings which are functions of g. Suppose 

f^\gl = ^m = -P{9^9*) (2.13) 

The above relation implies that if g{r) is a solution of the RG equations, then so is g*{rQ/r) 
for some tq. Furthermore, if g is self-dual at some scale tq, i.e. g = g*, then the ultra-violet 
(UV) and infra-red (IR) values of g are related by duality: giR = g^jy. 
The (3 function (|2.12| ) remarkably has such a duality with 

4 4 4r7 

9+ = —, 9- = —, 9a = 2- 2.14 

9+ 9- 9- 

The self-dual points g = g* are ((?+, g^, ga) = (±2, ±2, 0). 

The above duality is one ingredient we will use to resolve some difficulties encountered 
in ^ in completely interpreting the flows, as was done for anisotropic su{2) in [Q. But as 
we shall see not all RG flows may be resolved in this way. Also, unlike the su{2) case, we 
have not found an RG invariant for the above (3 functions. 



C. Virasoro master equation. 

The action ( ^.9| ) has a well-defined classical stress tensor Tl^]^^^. Introducing a 2d metric 
g'^'^, then T^Jf*^*^ = —4n6S/6g'^'^. One finds that this classical stress tensor is traceless, i.e. 
g^'^T^^^ = and the left-moving conformal stress tensor T = T^z is 

T:'Z'-9aT^ (2.15) 



iclass ^class 

' wzw 



where the T^ are defined in eq. ( |2.10|) , and T^^^ is the classical stress tensor for the WZW 



model which is the fc — >■ oo limit of the affine Sugawara stress tensor. 

In the quantum theory there are order g"^ corrections to the trace of the stress tensor and 
the conformal invariance is thus broken. Based on the form ( p.l5|) it is natural to suppose 



that if the theory fiows under RG to a conformally invariant fixed point, then at the fixed 
point the stress-tensor takes the form 

T&^cd point = Tg^ - ST, with 6T = HaT"^ (2-16) 

where Tg^ is the affine Sugawara stress-tensor and Ha are constants which are unknown 
functions of the fixed point values of qa- 

The condition that a stress tensor T satisfies the Virasoro algebra, i.e. 

T(^)T(0)~^ + 4t(0) + ...., (2.17) 

z* z^ 



is the so-called Virasoro master equation studied extensively in |Tl[. It is known that 
solutions come in pairs (K-conjugation), such that the sum of the stress tensors for a given 
pair is the affine- Sugawara stress tensor Tg^,. Thus bT (and consequently Tfixed point) must 
be a solution of the master equation, which reads 

2kD^^hAhB + C^^hAhB + 2C^^hAhB = he (2. 18) 

where the coefficients D^^, Cq^ and Cq^ are the same as those determining the /? function 
defined in eqs.( p.ll| ). Though it is intriguing that the same coefficients C,C,D determine 
both the /3 function and the Virasoro master equation, we cannot show in complete generality 
that solutions of the master equation correspond to zeros of the /3 function. 

Using the C,C,D computed in Ref. 0] one finds 2 pairs of solutions to eq. ( p. 181) . 



The fact that the solutions have a simple interpretation serves as a check of the coefficients 
C, C, D. The first pair of solutions is 

Tosp{2N\2N), = \{t-- T+) , To = (2.19) 

Here Tosp{2N\2N)-i is the affine-Sugawara stress tensor for the full current algebra and To for 
the empty coset osp(2A^|2A^)i/osp(2A^|2A^)i. The other pair of solutions is 



T 



gl{N\N) 



1= l^iT" -T j , Tosp{2N\2N)i/gl{N\N)^ =T - -T° - -T+ (2.20) 



corresponding to the sub-current algebra gl{N\N)i of osp{2N\2N)i and the indicated coset. 
(Tg/1-i = Tg — T-j-i). The two terms in Tg;(Ar|7v)i correspond to the two independent quadratic 
Casimirsof c//(iV|iV) 0. 

One can form a level-1 representation of both osp{2N\2N) and gl{N\N) with the same 
field content, thus both Tosp{2N\2N)i and Tg;(jv|jv)i are equivalent to the free-field stress tensor 
Tfree of the ip , P fields. This leads to the identity 

Tfree = Tg;(7v|Ar)i = Tosp{2N\2N)i, =^ T — -T" — -T = (2.21) 



This is equivalent to the statement Tosp{2N\2N)i/gi{N\N)i = 0. 

It is important to note that the purely random u{l) with only Qa not equal to zero does 
not appear as a solution of the master equation and is thus not a consistent stress-energy 
tensor. In terms of currents, 

11 ^^ 

TuiDo = 2^" = 2^'' ^ ^ ^(-)Wif- = Hp-H^ (2.22) 

a=l 

where 

N N 

Hp=Y. n^-^ H^ = ll r+r- (2.23) 

These currents satisfy two m(1) subalgebras: 

Hf,{z)Hp{0) ~ -/^v-W^^lO) ~ -^ (2.24) 

Note that H{z)H{0) ~ reg., i.e. since H satisfies a m(1) current algebra at level-0, and T„(i)g 
is formally the Sugawara form for u{l)o. The level being zero is what prevents Tu(i)„ from 
being a good stress-tensor. 

Consider the class of stress-tensors 

Th = Tfree " hTu(i)o = -^^^ + (1 " ^)^«(i)o (2-25) 

The stress-tensor for the pure random u(l) is of this form with h = 2ga, and corresponds to 
the only zero of the (3 function. This stress tensor satisfies: 

2 / \ 2 

Th{z)Th{0) ^ ^[Th- hT^^i)^) , Tft(z)T„(i)„(0) ~ — T„(i)o(0), T„(i),(2;)T^i),(0) ~ reg. 

(2.26) 

The Virasoro algebra is spoiled by the h term in the first equation. We should thus include 
this c = generalization of Virasoro as a possible fixed point. 

As far as the computation of conformal scaling dimensions is concerned we can formally 
view the stress tensors ( p.25[ ) as the coset '' osp{2N\2N)i/u{1)q\ However it is not a true 
coset since the m(1)o does not commute with Th. The point /i = 1 is special, in that Tu(i)q 
is precisely canceled and Th=i reduces to — T^/2, and this leads to a constant density of 
states. (See section IV.) Since the notation is convenient and descriptive we will refer to the 
Th=i theory as osp{2N\2N)i/u{1)q. We emphasize that the latter theory is not the same as 



the c = —2 theories osp(2N\2N)i/ u{l) ®u{l) considered in p!3|jr^JT0| J^ where the u(l)'s are 
generated by H^ and if/j. 

D. Singular RG trajectories in the Chalker-Coddington net^vork model. 

There exists a small category of RG trajectories which cannot be extended consistently 
to arbitrarily large RG time. This is probably linked to the level zero current algebras 
appearing in these disordered systems. 
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These pathologies concern trajectories, or their duals, which enter the perturbative phys- 
ical domain, the latter defined by requiring that the variances g^, gm, ga are all positive. In 
a finite RG time these singular trajectories converge to the singular locii g± = 2 with ga > 0. 

Let us first argue that the perturbative physical domain is stable under renormalization 
until reaching the singularities at g± = ±2. The boundaries of this domain are ga = 0, g-^- = 2 



and g+ ± g- = 0. From eq. (|2.12| ), one verifies that f3g^ > for ga > and < \g±\ < 2. 
Hence, ga remains positive as long as the RG trajectories are in this domain. Similarly, one 
verifies that the scalar product jSg-n, with n the normals to the boundaries pointing inwards, 
is positive on the boundaries. Thus, RG fiows which are initially inside the domain cannot 
cross the boundaries and stay within it. 

Since g^ increases, the only direction in which trajectories entering the physical domain 
from elsewhere in the phase space could go is towards the singular point g± = 2. Further- 
more, since in this domain ga is always increasing these physical trajectories reach g± = 2 
with ga > 0. This typical behavior may be checked by looking at the asymptotic behavior 
in the vicinity of the pole {g^,g-) = (2 + e_|_, 2 + e_). Here the dominant terms in the 
/3 functions are 

ft,.-% A_^5, ,3,,^% (2.27) 

Integrating these one finds 

ga^C exp e_/2, e+ ~ C exp 2C/e_ (2.28) 

with C and C two integration constants. Eq.( |2.28| ) shows that ga fiows to a non- vanishing 
constant as e_ -^ 0, while e+ vanishes exponentially. Since eq.( p.28|) involves two free 
parameters these trajectories are generic enough in this domain. 

Unfortunately, because these physical fiows reach the singularity at points distinct from 
the self-dual point {g± = 2,ga = 0), they cannot be extended smoothly beyond this point 
using duality and possible additional identifications as explained in p. Thus these RG 
trajectories are singular at a finite RG time, ie. after a finite scale transformation. This 
pathological behavior is in contradiction to what is expected for a renormalizable field theory. 
The singular fiows are depicted in figure 3. Trajectories can enter the physical domain either 
from the g^ = plane or from the point g± = ±2, ga = 0. 

E. Regular RG trajectories in the Chalker-Coddington network model and phase 

diagram. 

Most of the fiows are actually not singular. As was done for the anisotropic 5^(2) model 
in ref. |^, one may use duality arguments to construct smooth extensions such that these 
flows may be extended to arbitrarily large scale. The main steps in resolving these flows are: 
(i) RG flows to the singular self-dual points are extended beyond it using the duality and, 
(ii) flows where the couplings blow up in flnite RG time but do not cross self-dual points are 
extended with certain topological constructions with the (3 function serving as a consistent 
patching condition. In this section we describe the typical behaviors of these trajectories 
and identify the resulting asymptotic phases. As we will show, for this class of RG flows, 
there are only two such directions referred to as dXY and gX below. 

9 



To begin describing these regular flows, consider the (3 function when (7+ = 0. There, 

^^ = 0, Pg^ = -Agl/{2 - g_f (2.29) 

Since (3g_ = 0, we can consider the flows as originating or terminating at (7+ = 0. On the 
g+ = Q plane, the flows all originate from Qa = +C)0 since there Qa always decreases. Near 
(7+ = 0, one has 



l %+(5'«- 9ao) 
^2-9- 



^.^ 7:''':r . ^«o-^-(^--2)/4 (2.30) 



Thus, if Qa > Oao the flows are away from the g^ = plane, otherwise they are attracted to 
it. Since f3g^ < and g^ decreases in this domain, trajectories which are attracted to the 
plane (?+ = are more and more attracted to it, while those which flow away from it may 
flow back toward it if ga — gao is not large enough. This implies that some trajectories form 
loops, in the sense that they begin and end at different points on the g+ = plane but with 
different values of (?_ and ga- (See Figure 1.) 

Near g^ = oo one flnds Pg_ = 0, Pg^ = — 4(?^/(2 — (?_)^ and 



i2-g^y 



p,,^- '::'\ 'r\ gaoo^ {9^-2)12 (2.31) 



Thus if ga > gaoo, 9+ decreases away from oo, otherwise it continues to flow off to inflnity. 
Since Pg^ < 0, ga decreases there and trajectories flowing off will escape more and more to 
infinity, while those with g+ decreasing may be repelled back to infinity if ga — gaoo is not 
large enough, again forming loops that are simply duals of the loops near (7+ = 0. (See 
Figure 1.) Notice that gao{9-) = 9aoo{.9*-) in agreement with duality. 

Flows that do not loop back to (7+ = 0, 00 eventually encounter the poles in the /3 function 
in a finite RG time. Let us consider the ratios dg+/dg- = Pg^/Pg_ and dg+/dga = Pg+/Pga- 
They determine the direction of the fiows. One finds that, except at the self dual points, 
both dg+/dg^ = and dg+/dga = when g+ approaches the singularities at (?+ = ±2. This 
implies that the fiows near the poles at g+ = ±2 are tangent to the g^ = ±2 planes since 
P+ -C P- ,Pa- Similarly dg_/dg+ = and dg^/dga = when g_ = ±2. Thus /9_ <^ /9+ ,/?„ 
and the fiows are tangent to the g- = ±2 planes near g^ = ±2. Consequently if the fiows 
indeed cross the singular planes, they must do so at the self-dual points of g±. In order for 
the fiows to continue smoothly through the poles one also needs dg^/dg- and dgj^/dga to 
be finite at the poles. Examining the latter one finds that this can occur only when ga = 0. 
We thus conclude that if the fiows indeed cross the singular planes at finite angles, they 
must do so at the self-dual points g = g* for all g simultaneously. 

The crossing of the singularities may be described more precisely. Consider first the pole 
at g± = ±2. Looking for fiows approaching the self-dual point g± = ±2, ga = 0, one may 
search for an asymptotic expansion oi 2 + g^ and ga as functions of e+ = (?+ — 2. To lowest 
order, the equations for the trajectories are: 

g- = -2-e+ + 2Q_eX-Q_el + Oiel) (2.32) 

ga = Q-el-Q-el/2 + 0{el) 
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with Q_ a free parameter. These flows form a continuous set of trajectories. However since 
this is only a one parameter family, they are not the most generic ones. The (3 function for 
e+ is smooth, /3g^ = — 2 — e+ + • • ■ implying that e+ decreases. So these flows, coming from 
the domain g_ < 2, g^ > 2, reach the self-dual point in a flnite RG time and then cross 
over. (See Figure 2). 

The flows which cross the upper self-dual point ai g± = 2, Qa = either start with g+ < 2 
or g^ > 2. In the flrst case, their trajectories are described by 

^_ = 2 + e+ + Q+el + Q^e^ + 0{el) (2.33) 

ga = Q+el/2 + Qlel/2 + 0{eX) 

with 6+ = (7+ — 2 and Q+ a free parameter labeling them. The /? function Pg^ possesses a 
double pole, /3g+ = 32(1 + (Q+ - 2)e+ + • • •)/£+ so that e^ = 96(r - r*) + ■ ■ -. These flows 
reach the self-dual point from (7+ < 2 in a finite RG time but admit an extension beyond 
that point. As depicted in Figure 3 they then flow off to inflnity. In the second case, the 
trajectory equations are 

^- = 2 - e+ + (4g+ + 1)6^2 + • • • (2.34) 

9a = Q+e+ + ■ ■ ■ 

with g^ = 2 + e+. Again the /3 function for e+ has a double pole Pg^ = —32(1 + (4Q+ + 
l)e+ + ■ ■ ■)/^+ but with opposite sign. Hence, e+ now decreases and, after going through the 
self-dual point, these trajectories flow towards the axis g+ = 0. (See Figure 2.) 

Let us now show how one may resolve the blow-ups in flnite RG time. This concerns 
flows toward g-,ga = 00 that are dual to those crossing the g+ axis. Indeed, any trajectory 
that passes through (?_ = is mapped by duality to two separate trajectories at g- = ±cxd. 
Their asymptotic at inflnity are 



2 \2 



g+ -^ const., ga ~ const'. 5-!., g_ ~ -^- — ±^ (2.35) 

8gi{n - r) 



where r* is a constant. Eq.( ^.35D shows that in a flnite time r^, the flow is to {g^,g-,ga) 



(const., 00, 00). We have to extend the flows to larger scales. Since the trajectories are 
continuous at g^ = 0, it is natural to demand that their dual be continuous at (?_ = ±00 
so that flows towards {g-,ga) = (00,00) continue at {g-,ga) = (—00,00). Imposing this 
continuity requires identifying points at inflnity by changing the topology of the phase space 
in a way consistent with the (3 functions. To be more precise let us map points at inflnity to 
flnite distance by deflning new coordinates h±, ha by /i+ = g+, h^ = 1/g- and ha = —ga/d-- 
We identify the points with coordinates (/i+, ha, h^ = 0+) and (h^, ha, h^ = 0^) with ha 7^ 0, 
or alternatively 

{9+, 9-, 9a) = (9+, -9-, 9a), as \g^\ -^ 00, ga/gl ^ flxed (2.36) 

This is consistent with RG flows as the /3 functions Ph = Pg (dh/dg) in the new coordinates 
are continuous when /i_ crosses the origin. This proposal was described in greater detail for 
the su{2) case in |^. The identiflcation ( |2.36| ) changes the topology of the space of couphngs. 
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We can view this as gluing multiple patches of coupling constant space along g^ = ±00, and 
we will draw the resulting diagrams with this in mind. 

Having determined how regular flows navigate around singularities, we can list the tra- 
jectories which can consistently be extended to infinite RG time. These have been checked 
numerically, implying that they are generic enough. The flows are depicted in Figures 1-3. 
We classify them according to where they start from and where they terminate. The initial 
and final locations will either he g+ = or g+ = 00 as it takes an infinite time for the flows 
to reach these locii. For each flow there is a dual one, and we only list one of the two flows 
in the dual pairs. It is possible to draw the flows unambiguously in the g^, g^ plane since 
ga flows in a simple manner: in the UV ga = +00 and in the IR ga = — cx). 

Flows originating from the (7+ = plane: 

(a-1) The trajectory returns to g+ = at some other point without crossing the singularities 

and without g^ going to 00. (Figure 1.) 

(a-2) The flow is to g-,ga = cc with g+ < 2 without crossing the self-dual points. After 

identification of points at infinity as in eq. (|2.36|) , the flow then continues at g^ = — cxo and 

again eventually returns to the (7+ = plane. (Figure 2.) 

(a-3) The flow begins with < g- < 2 and first flows towards the self-dual point g± = 2, 

ga = 0. It goes through this point according to eq.( ^.34| ). After crossing the trajectory is 



again determined by duality, and flows to g^ = 00. (Figure 3.) 

Flows originating from the (7+ = 00 plane, dual to those terminating at g^ = 0: 

(b-1) The flow starts at (?+ = 00 and loops back to (?+ = 00 without crossing singularities. 

(Figures 1,3) The flow of this kind in Figure 3 is the dual of (a-2), and those in figure 1 the 

dual of (a-1). 

(b-2) The flow goes to g-,ga = 00 without crossing singularities. After the identification 

( p.36| ) the flow goes through the g± = ±2 poles according to eq. ( |2.32| ) and then to the 

g+ = plane. (Figure 2.) 

(b-3) The flow starts at infinity with < (?_ < 2, goes up to the self-dual point g±. = 2, 

ga = 0, crosses it according to (|2.34|) and terminates at g+ = 0. (Figure 2.) 



(b-4) The flow starts at infinity with g_ < —2, goes through the lower self-dual point at 
g± = ±2, ga = 0, according to ( |2.32|) , and then continues up to g+ = 0. It is a self-dual 
trajectory. (Figure 2.) 

All of the above flows asymptotically reach either 5f_|_ = or (yf_|_ = 00 at large scales. 
They take an infinite time to arrive at these locii. We shall refer to these phases as dXY 
and gX phases, respectively. The flows for (?+ < are the reflection of the flows for g^ > 
about the g^ axis due to the symmetry of the (3 function. 

dXY phase: g^ -^ 0, g- flows to a non-universal constant, and ga flows to —00 loga- 
rithmically, ga ~ — alogr with a = Ag'^/{2 — g-Y- Since the non-zero couplings {g^,ga) in 
the deep IR couple only the currents of the gl{N\N)i sub-current algebra this phase is in 
the same universality class as the disordered XY model, or Gade-Wegner class, as described 

rUf . Due to the structure of the P functions (|2.29|) since g- does not flow to 00, one 



m 



cannot argue that this phase flows to the fixed point osp{2N\2N)i/ gl[N\N)i as the solution 
of the Virasoro master equation would suggest, and it remains unclear whether one can 
identify a true IR fixed point to this flow. We point out that it was shown in |T^ that the 
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c = —2 coset conformal field theory osp(2A^|2A^)i/-u(l) ® u{l) describes a conformal sector 
but precise arguments identifying tliis witli tlie fixed point are missing. 

gX phase: Here, (/_ flows to a constant, Qa flows logarithmically to — oo, g^ ~ — alogr 
with a = 4(7^/(2 — g_Y whereas g^ grows much faster as (?+ ~ exp (2(7^/(7^). Since all the 
osp(2A^|2A^) currents are involved in this flow, this is probably a massive phase. 
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FIG. 1. RG trajectories that loop back to g^^. = 0,oo and flow to the dXY and gX phases 
respectively. 
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FIG. 2. RG trajectories that cross singularities and flow to the dXY phase. 
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FIG. 3. RG trajectories that either cross singularities and flow to gX phase or are singular. 
The physical domain is the 90° cone opening to the right. (Singular flows shown in red.) 



III. THE SPIN QUANTUM HALL NETWORK MODEL 

The hainiltonian of the network model for the spin quantum Hall transition is 

_ . 2a- a + m -idj + A- 
-idz + Az 2d ■ a — m 



where d,m are random potentials and A are random su{2) gauge potentials 
conventions for the Pauli matrices are 



a 



1/1 

2U -ly 



a 



J_ /O 1 

v^vo oy 



a 



J_fOO 
V2\l 



(3.1) 
Our 

(3.2) 
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With the above convention, Tra'^a'' = 5"^ /2 and /«'"^/«*'^ = -25^'^ where [a"", a^] = f'^'a^. 
The random potentials, m, a, A are taken to have the centered gaussian distributions with 
variance gm, ga, 9s- 

A. Effective action and its symmetries 

Introducing bosonic ghost partners P±,P± to the fermions, the effective action upon 
disorder averaging is 



where 



72 

SoS = 5frce + / ^ idaO'^ + gr^O^ + gsO-') (3.3) 



2 



O™ = i (^_V^+ - ^_^+ + /5_^+ - :5_/5+)' (3.4) 

In order to more clearly display the sub-algebraic structure, let us define new couplings 

(7,0" + (7„C>"^ = (7eO^ + (780^ ga = gc + ^gs, gm = \gs-gc (3.5) 

where (9^ = O^ — O"^, and O^ = ^0°' + ^O"^. The physical domain is the intersection of 
gc > —(78/2 and gc < 3(78/2 with gs > 0. The free conformal field theory has a maximal 
osp{4\4)k=i current algebra symmetry. The Lie superalgebra osp{4\4) has two commuting 
subalgebras osp{2\2) and su{2). The above notation reflects the fact that gs (gc) couples the 
spin (charge) degrees of freedom. Let us describe the osp{4\4) currents in a notation that 
clarifies these subalgebras. Define bosonic and fermionic su{2) currents 

L; = ^_a>+, L^ = p.a'^P+ (3.6) 

These satisfy level k = 1 and level k = —1 current algebra respectively 

^;w^Uo)~-^ + -r^^L^(o) (3.7) 



^2 2 z 

1 6"^ 1 



n{z)L'M -— + -r^^L^(o) (3.^ 



The operator O'^ is the su{2) invariant built from the sum: 

0^ = ^£"r, C = L^ + Lp (3.9) 

a 

The current C"" satisfies a level A; = current algebra. 

The currents in O^ generate an osp(2|2) subalgebra. Define 
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i id 

i i i,j 

where e^j = —eji, eu = 1, and similarly for the right-movers J = ip_^ip^,H = /?_,_/?_ 



These currents generate an osp{2\2)k current algebra at level k = —2 |2^. The currents 
J and J± generate a charged su{2) subalgebra under which the fermions ipi: transform as 
doublets. 

The operator O"^ is the osp{2\2) invariant bilinear: 

O^ = -Jl + HH + UJ-J+ + J+'J~) + S^^+ - S+S. + S+'S- - SJS+ (3. 10) 

To show this we have used the identity ( pTT|) of the Appendix. 

The currents in O^ are orthogonal to the above ones. To better reveal the su{2) properties 
of these currents, let us introduce intertwiners between the symmetric tensor product of two 
spin 1/2 representations and the spin 1 representation. These are: 

p- = a% ff = ea^, '={-1 o) ^^'^^^ 

We therefore define 

Bl = ^/5+P«/5+, 5- = ^/^-P'^/^- (3-12) 

Ul = (3+p^'^+, f/! = /?_p>_ (3.13) 

V^ = ^.a''(3+, V^ = (3_a''tlj+ (3.14) 

They transform in the spin 1 representation of su{2). The operator O^ can then be written 

as 

o^ = -2 (bJb+ + 5+5_ + uJj+ - f/+I7_ + y_F+ - y+F_ + LpTf^ - l^l^) (3.15) 

where e.g. BJB+ = Y.aB'^jB^. 

We can now describe the osp{2\2)®su{2) transformation properties of the above currents. 
The osp{2\2) currents ( |3.10|) commute with the su{2) currents C^. The currents in O^ 
transform in the spin 1 representation of su{2). The fields L^^p are not primary however: 

£'^(z)4(0)~l^ + -r^L^ (3.16) 

C\z)Ll{S))^--— + -r'^Ll 

Under the osp{2\2) symmetry the currents in O^ transform in the 8 dimensional inde- 
composable representation. This can be expressed as the OPE 

5'^(z)jr(0) ~ ^(t")}^'" (3.17) 

where S*" are the osp{2\2) currents, Jg'", i = 1, .., 8, a = 1..3 are the 24 currents in O^ with 
a the su{2) index, and t" is an 8 dimensional matrix representation of osp{2\2). 



IV. THE /3ETA FUNCTION FOR THE SPIN NETWORK MODEL 

The interactions that perturb away from the conformal field theory are of the form ( p.9|) 
where J" are the osp(4|4) currents at level k = 1 described above, and cif^ are quadratic 
forms that define the operators O"^ = {O'^, O'^, O^}. As before, the /3eta function proposed 
in ||2y] is expressed in terms of the data Cq^,D^^,Cq^ which can be computed from the 



OPE's (13.91). One finds the non-zero values 



Cf = cf = 4, Cf = 

yc8 ri8c 



Cf = cf = cf = -8 (4.1) 

The coefficients Cf and Cf are easily understood as minus the quadratic Casimir in the 
adjoint representation of su{2) and osp(2\2) respectively. 

The other RG data involves OPE's of the purely left-moving operators ( |^.10[ ). One can 



easily disentangle the D from the C term by noting that the D term only arises from central 
terms proportional to k in the OPE. One finds 

Df = -2, Df = Df = Df = 1 (4.2) 

Some of the C coefficients follow from simple group theory. Consider first OPE's with T^. 
Since T* has the form of the quadratic Casimir for su{2), we have T'^O^ ~ 20*, T^O^ ~ 20^ 
since 2 is the quadratic Casimir for the spin 1 representation of su{2). Thus Cf = Cf = 2. 
Similarly T'^O'^ ~ —40'^ where —4 is the Casimir for the adjoint representation of osp(2|2). 
The only non-zero OPE's of T"^ with O^ currents are 

T%z)L^ ~ -T%z)Lp ~ ^£ (4.3) 

which implies Cf = 8. The remaining C are computed tediously but straightforwardly 
using the free-field OPE's ( |2.6|) . For convenience we collect all non-zero C 



Cf = cf = 2 




f<88 _ _(^8s _ 


/-ICC 

^c 


Cf = 8, Cf 


= -6 



-c 



88 



(4.4) 



Using the general expression in |20|, the result for the beta functions is 



^- = (,, + 2)3(1 - 2) (^«' + ^«^^ - '^« + '^' 

^3= = (^1 _ 4)2^(^^ _ 1)2 {^^aWc - atal - i^ghl - Wc + 12^|) (4.5) 

^- - (,.-i)(,3lV(.8 + 2)^ (^-^^- - '^^'''^ + ''^^ - '''^ + '''' - ''^ 

+ ^9l9s{g8 - 2)igs + 2f{g, - 1) - 2g^{g^ - 2){gl + 4){g^ + 2f{g^ + 2g, 
The (3 function remarkably also has the duality ( |2.13|) with 
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9c = -, % = - 9s = ~-T (4-6) 

9c 98 98 

The self-dual points are {9s, 9c, 98) = (0, ±1,±2). As for the u(l) network model, we have 
not found any RG invariants, unlike the su{2) case. 



A. Virasoro master equation and critical exponents 

As for the u(l) network model, let us describe the solutions of the Virasoro master 
equation ( |2.18| ). We obtain four pairs of solutions. First one is the trivial one: 



To = and Tosp(4\4)i = jT" — -T 

Two other pairs of solutions are given by: 

1 3 1 

Tospi2\2)_2 = "g-^"" ^^^ ^osp(4|4)i/osp(2|2)_2 = g^'' " -^T 
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Tsu{2)o = -T' and Tosp(4|4)i/sn(2)o = -^T' + -T" - -T* 
The last pair of solutions, 

T=-T' - -r and T' = --T' + -T" - -T® 
2 8 2 8 2 

corresponds to osp{2\2)_2 ® su{2)q and to the coset osp(4|4)i/osp(2|2)_2 ® su{2)q. 

Let us recall the spin-charge separation of the free conformal field theory of the bosons 
and fermions ip, j3. This free theory has osp(4|4)i current algebra symmetry with stress tensor 
^osp(4|4)i- It was shown in Ref. [^ that the latter stress tensor decomposes as follows: 

Tosp{i\A)i = Tosp{2\2).2 + ^s«(2)o (4-7) 

where Tosp{2\2).2 ^^^ ^s«(2)o are the afiine-Sugawara stress tensors for the osp(2|2)_2 and 
su{2)q current algebras. This lead to the identity T' = 0: 

gyc _ 42-8 _4^x' = (4.8) 



It was argued in Ref. |2^ that though the stress tensor decomposes as ( |4.7| ), the Hilbert 
space does not exactly factorize due to logarithmic pairs such as 0^,0^^. This means that 
the coset osp(4|4)i/s-u(2)o is not identical to the current algebra osp(2|2)_2 even though it 
has a global osp{2\2) symmetry and the scaling dimensions are correctly given by osp(2|2)_2. 
However recent work shows that certain 4-point correlation functions can be factorized 



A comparison with other approaches was studied in ||24|| . 

The interesting possible fixed points are thus the cosets osp(4|4)i/su(2)o and 
osp(4|4)i/osp(2|2)_2. For each fixed point, let us compute some of the critical exponents. 
The density of states is proportional to the one-point function of $£;, 



piE) oc i^e) (4.9) 
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where ^£ is defined in eq. ( p.3|) . Let Te denote the scahng dimension of $£• in the JR. Then 

-^ oc E^s/(2-rB)^ asE ^0 (4.10) 

Other exponents are generally related to other perturbations $5: 

5S = 5[Sx^s (4.11) 



where 5 = is the critical point. Letting F^ denote the scaling dimension of $5 in the IR, 
the mass dimension of 5 is then 2 — Ts- As 5 ^ 0, there is thus a diverging length ^, 

e«5-^ ^ = 77^ (4.12) 

For a coset fixed point QkfHk', the dimensions of the fields in the IR are given by 

T = 2A = 2 (A^" - A^"') (4.13) 



The spin-charge separation implies that the conformal scaling dimensions A satisfy the 
sum rule: 

^osp(4|4)i ^ ^osp(2|2)_2 _j_ ^su{2)o (4.14) 

The osp(4|4) scaling dimensions are simply given by the free theory and are thus integer or 
half-integer, where ip,P have A^'^^*^'^'^-'^ = 1/2. The primary fields of the osp{2\2)_2 can be 
labeled (j, 6) where j is the spin of the su{2) generated by J, J^ and b is the m(1) charge 



H/2. These fields form 8j-dimensional multiplets and have left/right scaling dimension |25 



,2 ^2 



^oMm-. ^ J___ (4^^5) 

The primary fields of su{2)o on the other hand are labeled by a spin j and have scaling 
dimension [^ 



osp{4\4)i/osp{2\2)_2 The osp(2|2)_2 dimension of $£■ is 1/4 since it transforms in the 
(1/2, 0) representation with A""p(2|2) ^ ^/g_ r^^^^g r^ = 1 - 1/4 = 3/4 and one finds 

]iE)(xE^/^ (4.17) 

In this phase it is the charge degrees of freedom that are massive, i.e. localized. 

osp(4|4)i/su(2)o. The field ^e has j = 1/2 quantum numbers under the su{2)q, and 
thus has A^"(2) ^ 3/8_ jj^jg ^g^^^g to F^ = 1 - 3/4 = 1/4 and 

pC^ocE^/^ (4.18) 

Let us compare the latter coset with percolation. The map to percolation involves 
averaging over su{2) matrices with the fiat Haar measure, which is expected to correspond 
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to gs = oo. The gs = oo theory was described in ||2^. This averaging projects out the su(2) 
degrees of freedom and in our approach this is analogous to arguing that the RG flow gaps 
out the sm(2)o sub-current algebra. We do not have an argument leading to the identification 
of $5 for the localization length exponent. However there does exist a field in our theory 
with the correct dimension to lead to the percolation exponent. Consider a field $5 in the 
free theory of the form $5 = <f)s4'5 with (ps ~ [tpipip (3 (3) . In the free theory this field has 
A = 5/2. Since ip,P have j = 1/2 under the su{2)q, we can take (ps to be in the j = 3/2 
representation of su{2)q with A''^^^^^ = 15/8. Under the osp{2\2)_2 this field transforms in 
the (j, b) = (3/2, 1) representation with A^^^^^'^^ = 5/8. Note that the sum rule is satisfied: 
15/8 + 5/8 = 5/2. Since Ts = 5/4 this leads to z/ = 4/3. 



B. RG flows 

In this section we describe the RG flows based on the (3 function ( [4.5| ). As in the u{l) 
network model, there exist both well-behaved and singular flows. 

General properties. 

When gs = the spin and charge degrees of freedom decouple and the P function reduces 
to 

The fact that the one-loop result for Pg^ is exact is due to the level k being zero for the su{2) 
currents C^. If gs > 0, gs flows to infinity after a finite scale transformation, whereas ii gs < 
it flows to zero, i.e. gs > is marginally relevant and gs < is marginally irrelevant. On 
the other hand ^fc < is marginally relevant with gc flowing to —00 and gc > is marginally 
irrelevant with gc flowing to zero. 

When gs 7^ 0, we can consider the flows as originating from the gs — gc plane. On 
this plane the flows originate in the UV from gs = 0, —00 or gc = 0, +00. Near gs = 0, 
/?gg = 2gsgs- Thus ii gs > {gs < 0) one flows away (toward) the gs — gc plane. If gs is 
not large enough, the flows can turn around and form a loop returning to the gs = plane. 
Otherwise it continues away from the plane until it reaches a singularity. By duality there 
are similar properties at f^g = cxd. 

At generic points, the slope of the trajectories dgs/dgs = Pgs/Pgs is zero at gc = ±2, 
whereas the slope dgs/dgc = Pgs/Pgc is zero at f^g = 2 but equals cxd at (^g = ~2. Thus 
the flows are generically tangent to the gs = 2 plane, since Pg^, Pg^ ^ Pg^ near ^fg = 2. 
Furthermore, Pg^ ~ 6(5'^ — l)^/((y'g — 2)^, so that gc increases in the vicinity of the plane 
gs = 2. On the other hand, since Pg^ -C Pgg <C Pg^ near gs = —2 on the other hand, the 
flows are generically attracted or repealed by the gs = —2 plane. More precisely, because of 
the odd order pole, Pgg ~ —?>2gs/{gs + 2)^, the flows are attracted to gs = —2 for gs > 0, 
whereas for f^s < flows are away from f^g = —2. One also has Pg^ ^ 16g^/{gs + 2)^ near 
gs = —2, the coupling gs very rapidly increases. 

The behavior near the poles f^c = 1 is somewhat simpler. There, the slopes are such that 
dgs/dgc and dgs/dgc both vanish since Pg^ -C Pg^ -C Pg^ at generic points. Near (7c = 1 and 
away from gs = ±2 the beta function Pg^ ~ —2/{gc — 1)^ has double pole. Thus gc always 
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decreases and flows can cross the Qc = ^ plane. Note that Qc = —1 is not a singular point. 
Near g^ = —1 one flnds (3g^ = —1/2, so again the flows easily cross the Qc = —1 plane. 

The non generic singular points are those for which the ratio of the beta functions is 
finite and non- vanishing. These are at {qq = 2,gc = ±1). At Qc = —IjQs = '^ and (?« = 
the ratio dgs/dgc and dgs/dgs are finite. Thus one can fiow through this self-dual point and 
extend the fiow by duality. We have numerically verified such fiows. As analyzed below, the 
fiows toward {gc,g%) = (1,2) on the other hand are singular. 

As in previous examples, blow-ups in finite RG time r are resolved by certain topological 
identifications. When g^ ^ 0, there are fiows that blow up in finite RG time in the direction 
(7c ^ oo. As can be seen numerically, in this domain g^ increases rapidly while g^ fiows close 
to the value g^ = 2. This can be seen by expanding near gs = 2 + es with eg — > 0. The 
/3-functions asymptotically reduce to 

A^-^, (3c ^^, /38^-^ (4.20) 

^0 ^8 ^8 

The solutions are 



gs ^ Ve8 + const., eg ~ 2c^(r* - r), 9c^ - ^ log((r* - r)6) (4.21) 

showing that (gs^gcds) ~^ (const., oo, 2). Note also that gc blows up before eg vanishes. 
To continue the fiows note that the asymptotic behavior ( [4.20| ) satisfies P{gs,gc,g8) = 
Pigs,~9c,g8) as \gc\ — > C)0. Viewing this as a consistent patching condition, we can thus 
identify gc = +oo with gc = — oo when gg ^ 0, such that fiows toward gc = oo continue at 
gc = -oo. 

Singular fiows. There are two classes of singular fiows: 

(i) Consider the fiows toward the self-dual point gc = l,g8 = 2. Letting gc = 1 + ec, 
gs = 2 + es, asjTiiptotically the /3 function behaves as 

^2 32 16-4^^-^2 24 

/?g^ ~ -^ + + ^ — ^ + _ + ... 

^8 ^c^S ^8 Cc 

24 2 12 4 

/3..-^-4 + --- + --- (4-22) 

4 ^c eg ec 

/?.s--— + ■•• (4.23) 

The trajectories fine tune themselves coherently in order to converge toward the singularity 
and to cancel the double poles of the /5-functions. The RG trajectories may be found as 
series expansion in terms of eg, with eg — *> 0. We obtain 

ec = eg(l + a^eg + ■ • •)/2V3 

with ttc solution of 3(4v^ac + v^ — l)/2 = —3^/'* and 

^, = 8- 3^/^(1 + ^-i^eg log eg + ---) 
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The RG equation for eg, eg = Ps, then gives 



eg ^ -3i/V32(r, - r) + ■ ■ ■ 

showing that the trajectories are effectively attracted by the pole. Since gs flows to 8 ■ 3^/^ 
rather than zero, this flow cannot be continued through the singularity by duality. These 
singular flows are in the physical domain. 

(ii) Consider the flows toward gs = —2. Letting g^ = —2 + eg, the /? function reduces to 

^8 ^8 ^8 

Integrating these equations one finds the asymptotic behavior 

9s^ces'^\ ^^A,f + c', es ^ {lUcin - r)f/' (4.25) 

gc + i- 40c 

for some constants c,c',r^. Thus in a finite RG time {gs,gc,gs) —^ (cxd, const., — 2). Again, 
there is neither the duality nor a consistent patching condition allows the flows to extend 
beyond this singularity. 

Regular flows All other flows can be extended to arbitrarily large or small length scales. 
These well-behaved flows either don't encounter singularities or can be extended consistently 
by passing between patches patches glued at gc = ±cxd, and/or flowing through the self-dual 
points {gs, gc, gs) = (0, —1, 2). The flows that do not loop back to (^g = or cxo are shown in 
Figures 4,5. 

Below we list the asymptotic directions, or phases, that can be continued to arbitrarily 
large length scales. 

phase O In this phase all disorder is marginally irrelevant. The flows are attracted to 
gs = with gs < 0,gc> and then gs, gc continue to flow to zero according to the (3 function 

(EH)- 

phase C This phase is similar to phase O except that gc < and thus marginally 
relevant. Thus gs,g?, ^ ^ and gc -^ — oo. 

phases S^ A large set of flows, some with rather different UV properties, asymptotically 
flow in the direction S~^. Some flows originate in the UV with gc > 0. Flows leaving the 
gc — gs plane with gc > flow to ^^c = cc, continue at gc = — oo, through the self-dual 
point {gs,gc,gs) = (0,-1,2) and then in the S^ direction. These flows originate from 
{gs, gc, gs)uv = (O'*', — oo, O"*"), and then by duality should flow to g*, i.e. in the IR we expect 
{gs, gc, g^liR = (— oo, 0, +oo). As shown below, the asymptotic behavior of S~^ supports this. 
Flows originating from {gs,gc,g%) = (+oo,0, +oo) are also attracted to 5"^ but do not pass 
through self-dual points. 

Assuming gsgc ^ 1, the /3 function asymptotically reduces to 

Ps ^ -Sgs/gs, (3c ^ 2gl, (3s ~ -I2gl/gl - 8gs 

The equation for gc decouples and gc — 1/r with r the RG time. Integrating the RG 
equations in the gs, gs plane gives the equations for the asymptotic trajectories, 

gl ~ 2gl log 5-8, 98 oc r^ log r, gcOcl/r 
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for large r. As a consistency check one verifies that g^ ^ g"^ and gsgc ^ 1- The phase S 
has the same properties as S~^ with gs — > — oo. 

phase SC^ These flows originate from the gc — gs plane with g^ < 0. In the deep UV 
one has {gs, gc, g»)uv = (0"'', 0~, 0+). Since the flow is through a self-dual point, one expects 
from duality that {gs,gc,g%)iR = (+oo, -cx), +oo). Assuming \gs\ < \gs\ and l^ffl < Ifi'cfl'sl, 
the /3-functions reduce to 

13% ^ , 13c ^ -2, [5s ~ — 

9s gi 



The trajectories read 

for large r where c is some constant. The phase SC~ is very similar to SC~^ but with gs < 0. 



9^s - cgl, gs ^ 16cr^ gc ~ -2r 
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FIG. 4. Flows away from gs = m the spin network model. 
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FIG. 5. Flows toward gs = in the spin network model. 

C. Identifying the fixed points 

It is possible to interpret the asymptotic directions listed in the previous subsection as 
some of the above coset solutions of the master equation. The main hypothesis in making 
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this identification is the following. Suppose a coupling g flows to infinity much faster than 
the other couplings and g couples a sub-current algebra Hk'- Then we argue that the Hk' 
currents are gapped out in the flow leading to the fixed point Qk/'Hu'- 

Phase C Since the only non-zero coupling is gc and gc -^ — oo, this phase clearly has 
the interpretation as the fixed point: 

Phase C IR fixed point : osp{4\4:)i 

Phases S^ . In the phases S^ the only couplings flowing to cx) are g8,gs- It was shown 
in Ref. |]22| that the operators O^, O'^ are a logarithmic pair. One can see explicitly from the 



asymptotics that gs/gs flows very rapidly to zero. We can thus interpret this as gg effectively 
being zero in comparison to gs- Since gs is flowing to ±00, we propose the fixed point is the 
coset osp(4|4)i/sm(2)o. 

Phases SC^ . In this phase gc decouples in the RG equations and flows to —00 with the 
same asymptotic as phase C. Hence both subalgebras su{2)o and osp(2|2)_2 are massive, 
and we propose the flxed point osp(4|4)i/su(2)o®osp(2|2)_2 = 0- Since this theory is empty, 
this is a massive phase. This leads to P^; = and a constant density of states as -E — > 0. 

V. MULTIFRACTAL EXPONENTS FROM CONFORMAL COSETS 

Critical exponents related to the multi-fractality of the wave functions have been the 



subject of numerous studies [p^p8| - |30|JT6[| . The multi-fractality itself is an important con- 



straint on the critical theory. These exponents and the correlation length exponent C, oc E 



where u ~ 2.35 ||3l]] are both important benchmarks for the correct critical theory. The 
multi-fractal exponents are more easily related to the conformal scaling dimensions of fleld 
operators in the critical theory [Q |I8| . 

In this section we show that the critical theories considered above lead to multifrac- 
tality in a very natural way. For the original m(1) network model we will consider the 
osp{2N\2N)i/u{l)o conformal fleld theory deflned in section IIC, with stress tensor 

T = Tosp{2N\2N)i — Tu{l)o (5.1) 

For the spin network model we will consider the coset osp{AN\AN)i/su{2)q. 

For any fleld $(x) in an effective fleld theory of disorder we can deflne the quantity 



P(^) = /dM^(-))l (5.2) 



where ($(x)) denotes the disorder average of the correlation function ($(x)), and {^{x))i 
denotes the disorder average of its q-th moment in the effective theory with N > q copies. 
The above scales as 

p(i) ~ L-^(i) (5.3) 
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where L is a large length scale. From simple dimensional analysis one finds 

r(g) = r(g)-gr(l) + 2(g-l) (5.4) 



where T{q) is the scaling dimension of ($)^. 

If r(g) is quadratic in q, then imposing r(0) = —2 and r(l) = 1 leads to a one-parameter 
family 

r(g) = (2 - «o)g' + «og - 2 (5.5) 

It is known that usually r(g) changes form for q above some threshold [^. However we are 
interested in the parameter oq which is determined near q = since ao = dr /dq\q=Q, so we 
do not consider the phenomenon of multi-fractality termination. It is customary to define 
the Legendre transform with a = dr/dq, which defines q{a), and 

/(a.) = a,-r(,) = -|^ + 2 (5.6) 

The fact that the numerical studies show a close fit to the above parabolic form indicates 
that r(g) is indeed quadratic in q for small q. The exponent ao determines the typical 
amplitude 

Atypical = exp (M$y) ~ ^~°° (5.7) 

The multi-fractal dimensions T{q) follow simply from the conformal dimensions with 
respect to the current algebras using eq. ( |4.13| ) where Q = osp{2N\2N)i and H = u{l)o or 



sm(2)o, and A^'^ are conformal scaling dimensions which follow from the operator product 
expansion with the stress tensor. 

It is convenient to bosonize the fields, where for each copy 

^± = e±^^ /9+ = e^^'r/, /3_ = e'^'^'^.e (5.8) 

where (p, ip' are the 2;-dependent parts of scalar fields = ip{z) +'ip{z), and 0' = f'{z) +Tp'(z) 
and (?7,0 is a conformal dimension (1,0) fermionic system with central charge c = —2 ||33|| . 
The bosonized action is then 

5free = y ^ ( E ^<f^^^'^^^ " Ocf^'M + V%C + tdX\ (5-9) 

The m(1) currents are 

N N 

H^ = J2^i^- = ^E'9.</>., Hp = T.P+P- = ^Y.dz(f)[ (5.10) 

i=l i i=l i 

and similarly for H. 

Consider first the fields $ = exp(m0) in the ■u(l) theory labeled by a parameter a. T{q) 
is the dimension of 

[$], = $^(x)$2(x)---$5(x) (5.11) 
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where $* is the i-th copy of the at least N = q copy effective theory. The operator [$]g has 
dimension qa"^ with respect to the osp{2N\2N)i since each copy has dimension a^. Since the 
H u{l) charge is ag, eq. ( p.l| ) imphes 

V{q)=a\{l-q), => ao = 2 + a\ (5.12) 

for $ = exp{ia(f)). 

Note that it is important for the denominator of the coset to be independent of N in 
order to obtain this resuh. This is equivalent to the statement that T{q) can be computed 
at arbitrary N > q and should give the same result independent of A^. In particular the 
c = — 2 theories osp{2N\2N)i/u{l) ^^(l) where the u(l)'s are generated by H^, Hp do not 
have multifractal properties due to the A^-dependence in eq. ( p. 241) . 

The multi-fractal exponents for the wavefunctions \l/(a;) correspond to choosing # as the 
density operator 

p = ifj^ip^ + tp-ip _^ = CO?, (p (5.13) 

The reason is that |\1/|^ oc Gy:ct{,x,x) oc (p). Since r(l) = 0, this implies a constant density of 
states at zero energy. Since p corresponds to a = 1, this naively gives Oq = 3. However one 
should perhaps bear in mind that the numerical simulations of the network model used pe- 
riodic boundary conditions on the wavefunctions in a cylindrical geometry. It is well-known 
that periodic boundary conditions on the cylinder are mapped to anti-periodic boundary 
conditions on the plane. The boundary conditions of free fermion fields are modified in the 
presence of spin, or twist, fields, which here are a± = exp{±i(f)/2). This can be seen from 
the OPE 

which shows that ip± picks up a phase —1 as 2; encircles the origin on the plane. Since 
the more relevant operator in the OPE of p with the twist fields is again a twist field, we 
suggest that ao may correspond to a = 1/2, which gives oq = 9/4. This is quite close to 
the published simulations on the network model |1THITB[, which report Oq = 2.26 ± .01 and 



2.260 ± .003. The small errors in the latter measurement seem to rule out 9/4. Simulations 
of other models believed to be in the same universality class give the slightly higher result 
ao = 2.28 ±.03. 

The point-contact conductance is another multi-fractal exponent that has been studied 



numerically for the network model p9| , |30[| . Define higher moments of the density-density 
correlation as 

G(9) = (p(r)p(O))'? ~ r-"«(^) (5.15) 

where p is the density operator. Let us assume that in the flow to the IR there is no mixing 
between the various copies. Simple dimensional analysis then gives 

rciq) = 2r(g) (5.16) 



where F is defined the same way as above, i.e. as the dimension of ($)'?. 
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It is possible to derive a general scaling relation between ao and the typical point contact 
conductance exponent Xt. Let us generally suppose that r(g) = Aq + Bq^ for some constants 
A, B. If the density of states is constant, this implies r(l) = 0, or A = —B. Then, for the 
point-contact conductance, eq. (|5.16| ) leads to an /(a) spectrum 



4Xi 
where Xt = 1A is the typical point-contact conductance 

exp (logG) ~ r~"^' 



{a - Xtf 
Jg{o-) = tt; (5.17) 



Comparing with eqs. ( |5.4|J5.5| ) we have ao = 2 + A, which implies the scaling relation: 

Xt = 2{ao-2) (5.18) 



We expect the above relation to hold whenever the density of states is constant |[T^. The 
above relation works reasonably well using the most recent values of Xf of .54, .57 reported 
in [^ and the numerical values of ao PH| , P!3| , P!B| . For oq = 9/4, one has Xt = 1/2. 



Let us turn next to the coset osp{2N\2N)i/su{2)q. The fermion fields transform as spin 
j = 1/2 doublets under the su{2). Consider again the su{2) invariant density operator 
$ = p = ipip, with dimension 1 in the osp{4\4)i theory. Since su{2) primary fields of spin j 
have dimension A*"^ = j{j + l)/2 when k = |2^, one finds r(l) = 1/4. Since [$]q has 
spin j = q/2, we thus find 

r{q)=q-q{q + 2)/4 (5.19) 

This agrees with ref. |]l8l at small q. Interestingly, using eqs. (p.4| ,^.5| , ^^) , the above T{q) 
also leads to ao = 9/4. For the point-contact conductance, here the density of states is not 
constant, i.e. r(l) ^ 0, so eq. (|5.18| ) does not apply. Rather, eq. (|5.16| ) leads to Xt = 1. 



If we adopt the viewpoint that the above computations are a positive indication, it is 
tempting to try and understand the correlation length exponent z/ for the quantum Hall 
transition. This is less straightforward since this exponent is not associated with the density 
operator. Rather, it is expected to be related to a perturbation of the action 6S = 6 J d?x $5 
where tuning 5 through zero corresponds to tuning the network model to its critical prob- 
ability. As described in section IVA, this leads to a diverging length ^ oc 5'" , where 
v = l/{2-Ts). 

The challenge is to determine $5 from first principles and this is beyond the scope of 
this paper. Even for the case of the spin quantum Hall transition the arguments leading 
to the identification of the operator $5 which gives the percolation exponent z/ = 4/3 in 
this approach are missing, though there does exist an operator of the right dimension in the 
coset osp(4|4)i/s'u(2)o [^. (See section IVA.) In the percolation picture, $5 is the two-hull 
operator. In the remainder of this section we just explore the kinds of exponents that can 
arise. 

Since the 1-copy theory has a non-critical density of states we must consider at least the 
2-copy theory. Since the density operator is coscj), it is natural to consider 
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where $^ /^ is a field in the i-th rj/^ copy. Each 77/^ copy describes dense polymers with a 
variety of non-trivial dimensions [^, and dividing by m(1)o does not affect these scaling 
dimensions. Interestingly, there have been some recent attempts at describing the Quantum 



Hall transition starting from polymers rather than percolation ||3^. See also ref. [0,0] . 

Dividing by u(l)o endows the combined expi0 factors with dimension 2. If we simply 
take ^lijc as the 1-leg operator with V = —3/16 and ^'i,. as a 0-leg, or twist, operator with 
r = —1/4, this gives Ts = 25/16. This leads to z/ = 16/7 ~ 2.29. This appears to be the 
possibility that is closest to the numerical value of 2.35 in this scheme @]. 



VI. CONCLUSIONS 

We have described how the RG fiows based on the all-orders f3 function for both the 
Chalker-Coddington and spin network models are attracted to a true singularity after a 
finite scale transformation. The nature of these singularities remains to be understood. The 
RG fiows in other domains are regular and can be continued to arbitrarily large length scales. 

Using an argument that is independent of the P function, we proposed that the fixed 
points of general current- current interactions correspond to solutions of the Virasoro master 
equation. For the network models these correspond to coset conformal field theories, the 
most promising being osp{2N\2N)i/u{l)o for the u{l) network and osp{4N\4N)i/ su{2)q 
for the spin network. Both of these theories lead to multifractal properties that we have 
analyzed. 

For the spin Quantum Hall transition, it remains an open question whether the 
osp{AN\4N)i/ su{2)q for A^ = 1 is equivalent to percolation. This can perhaps be answered 
by comparing with the recent partition functions in [^. It would also be interesting to 



perform numerical simulations of the multifractal properties of the spin network model and 
compare them with the coset prediction of ag = 9/4, since this result relies on averages of 
all higher moments and is thus not accessible from percolation. 
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VIII. PAULI MATRIX AND INTERTWINER IDENTITIES 

Using the relations [a", a^] = /'^''V^ {cr^^a''} = ^S'"'' and e^ = -1, (a'^)t = a" , and 
ea'^ = —cr'^e, one can establish the following identity 

(^ij(^nm = 2^imSjn " -^SijSnm (8.1) 

The intertwiners satisfy the following identities 
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rp' = ^/"'v^ - i^"^ 



a>* 






P"^ 



H 2 4 

2^ ~ 4 



^.2) 
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